A WOX homolog disrupted by a transposon led to the loss of spines and contributed to the domestication of lettuce.

The loss of spines is one of the most important domestication traits for lettuce (Lactuca sativa). However, the genetics and regulation of spine development in lettuce remain unclear. We examined the genetics of spines in lettuce using a segregating population derived from a cross between cultivated and wild lettuce (Lactuca serriola). A gene encoding WUSCHEL-related homeobox transcription factor, named as WOX-SPINE1 (WS1), was identified as the candidate gene controlling the spine development in lettuce, and its function on spines was verified. A CACTA transposon was found to be inserted into the first exon of the ws1 allele, knocking out its function and leading to the lack of spines in cultivated lettuce. All lettuce cultivars investigated have the nonfunctional ws1 gene, and a selection sweep was found at the WS1 locus, suggesting its important role in lettuce domestication. The expression levels of WS1 were associated with the density of spines among different accessions of wild lettuce. At least two independent loss-of-function mutations in the ws1 gene caused the loss of spines in wild lettuce. These findings provide new insights into the development of spines and facilitate the exploitation of wild genetic resources in future lettuce breeding programs.

Lskipk Lsatpase double mutants are necessary and sufficient for the compact plant architecture of butterhead lettuce.

Lettuce, an important leafy vegetable crop worldwide, has rich variations in plant architecture. Butterhead lettuce, a popular horticultural type, has a unique plant architecture with loose leafy heads. The genetic and molecular mechanisms for such a compact plant architecture remain unclear. In this study we constructed a segregating population through crossing a butterhead cultivar and a stem lettuce cultivar. Genetic analysis identified the LsKIPK gene, which encodes a kinase, as the candidate gene controlling butterhead plant architecture. The Lskipk gene in the butterhead parent had a nonsense mutation, leading to a partial predicted protein. CRISPR/Cas9 and complementation tests verified its functions in plant architecture. We showed that the loss of function of LsKIPK is necessary but not sufficient for the butterhead plant architecture. To identify additional genes required for butterhead lettuce, we crossed a butterhead cultivar and a crisphead cultivar, both with the mutated Lskipk gene. Genetic mapping identified a new gene encoding an ATPase contributing to butterhead plant architecture. Knockout and complementation tests showed that loss of function of LsATPase is also required for the development of butterhead plant architecture. The Lskipk Lsatpase double mutation could reduce leaf size and leaf angle, leading to butterhead plant architecture. Expression and cytology analysis indicated that the loss of function of LsKIPK and LsATPase contributed to butterhead plant architecture by regulating cell wall development, a regulatory mechanism different from that for crisphead. This study provides new gene resources and theory for the breeding of the crop ideotype.

PSEUDO-ETIOLATION IN LIGHT proteins reduce greening by binding GLK transcription factors.

Knocking out genes encoding proteins that downregulate the accumulation of pigments may lead to increases in crop quality and yield. PSEUDO-ETIOLATION IN LIGHT 1 (PEL1) downregulates the accumulation of carotenoids in carrot and chlorophyll in Arabidopsis and rice and may inhibit GOLDEN 2-LIKE (GLK) transcription factors. PEL1 belongs to a previously unstudied gene family found only in plants. We used CRISPR/Cas9 technology to knock out each member of the 4-member PEL gene family and both GLK genes in Arabidopsis. In pel mutants, chlorophyll levels were elevated in seedlings; after flowering, chloroplasts increased in size, and anthocyanin levels increased. Although the chlorophyll-deficient phenotype of glk1 glk2 was epistatic to pel1 pel2 pel3 pel4 in most of our experiments, glk1 glk2 was not epistatic to pel1 pel2 pel3 pel4 for the accumulation of anthocyanins in most of our experiments. The pel alleles attenuated growth, altered the accumulation of nutrients in seeds, disrupted an abscisic acid-inducible inhibition of seedling growth response that promotes drought tolerance, and affected the expression of genes associated with diverse biological functions, such as stress responses, cell wall metabolism hormone responses, signaling, growth, and the accumulation of phenylpropanoids and pigments. We found that PEL proteins specifically bind 6 transcription factors that influence the accumulation of anthocyanins, GLK2, and the carboxy termini of GLK1 and Arabidopsis thaliana myeloblastosis oncogene homolog 4 (AtMYB4). Our data indicate that the PEL proteins influence the accumulation of chlorophyll and many other processes, possibly by inhibiting GLK transcription factors and via other mechanisms, and that multiple mechanisms downregulate chlorophyll content.

Seed color in lettuce is determined by the LsTT2, LsCHS, and Ls2OGD genes from the flavonoid biosynthesis pathway.

KEY MESSAGE: The mutated LsTT2 and Ls2OGD genes are responsible for white seeds and yellow seeds in lettuce, respectively. Three LsCHS genes are involved in the biosynthesis of flavonoid in seed coats. Lettuce seeds have several different colors, including black, yellow, and white. The genetic mechanisms underlying color variations of lettuce seeds remain unknown. We used genome-wide association studies (GWAS) and map-based cloning approaches to clone genes controlling the color of lettuce seeds. LsTT2, which encodes an R2R3-MYB transcription factor and is homologous to the TT2 gene in Arabidopsis, was shown to be the causal gene for the variation of black and white seeds in lettuce. A point mutation leads to the lack of stop codon in the LsTT2 transcript, resulting in white seeds. Knockout of the LsTT2 gene converted black seeds to white seeds. The locus controlling yellow seeds was mapped to Chromosome 2. Knockout of two 2-oxoglutarate-dependent dioxygenases (2OGD) genes from the candidate region converted black seeds to yellow seeds, suggesting that these two 2OGD proteins catalyze the conversion of yellow metabolites to black metabolites. We also showed that three LsCHS genes from the candidate region are associated with flavonoid biosynthesis in seeds. Knockout mutants of the three LsCHS genes decreased color intensity. This study provides new insights into the regulation of flavonoid biosynthesis in plants.

Folate shapes plant root architecture by affecting auxin distribution.

Folate (vitamin B9) is important for plant root development, but the mechanism is largely unknown. Here we characterized a root defective mutant, folb2, in Arabidopsis, which has severe developmental defects in the primary root. The root apical meristem of the folb2 mutant is impaired, and adventitious roots are frequently found at the root-hypocotyl junction. Positional cloning revealed that a 61-bp deletion is present in the predicted junction region of the promoter and the 5' untranslated region of AtFolB2, a gene encoding a dihydroneopterin aldolase that functions in folate biosynthesis. This mutation leads to a significant reduction in the transcript level of AtFolB2. Liquid chromatography-mass spectrometry analysis showed that the contents of the selected folate compounds were decreased in folb2. Arabidopsis AtFolB2 knockdown lines phenocopy the folb2 mutant. On the other hand, the application of exogenous 5-formyltetrahydrofolic acid could rescue the root phenotype of folb2, indicating that the root phenotype is indeed related to the folate level. Further analysis revealed that folate could promote rootward auxin transport through auxin transporters and that folate may affect particular auxin/indole-3-acetic acid proteins and auxin response factors. Our findings provide new insights into the important role of folic acid in shaping root structure.

Watermelon domestication was shaped by stepwise selection and regulation of the metabolome.

Although crop domestication has greatly aided human civilization, the sequential domestication and regulation of most quality traits remain poorly understood. Here, we report the stepwise selection and regulation of major fruit quality traits that occurred during watermelon evolution. The levels of fruit cucurbitacins and flavonoids were negatively selected during speciation, whereas sugar and carotenoid contents were positively selected during domestication. Interestingly, fruit malic acid and citric acid showed the opposite selection trends during the improvement. We identified a novel gene cluster (CGC1, cucurbitacin gene cluster on chromosome 1) containing both regulatory and structural genes involved in cucurbitacin biosynthesis, which revealed a cascade of transcriptional regulation operating mechanisms. In the CGC1, an allele caused a single nucleotide change in ClERF1 binding sites (GCC-box) in the promoter of ClBh1, which resulted in reduced expression of ClBh1 and inhibition of cucurbitacin synthesis in cultivated watermelon. Functional analysis revealed that a rare insertion of 244 amino acids, which arose in C. amarus and became fixed in sweet watermelon, in ClOSC (oxidosqualene cyclase) was critical for the negative selection of cucurbitacins during watermelon evolution. This research provides an important resource for metabolomics-assisted breeding in watermelon and for exploring metabolic pathway regulation mechanisms.

Loss-of-function of SAWTOOTH 1 affects leaf dorsiventrality genes to promote leafy heads in lettuce.

The mechanisms underlying leafy heads in vegetables are poorly understood. Here, we cloned a quantitative trait locus (QTL) controlling leafy heads in lettuce (Lactuca sativa). The QTL encodes a transcription factor, SAWTOOTH 1 (LsSAW1), which has a BEL1-like homeodomain and is a homolog of Arabidopsis thaliana. A 1-bp deletion in Lssaw1 contributes to the development of leafy heads. Laser-capture microdissection and RNA-sequencing showed that LsSAW1 regulates leaf dorsiventrality and loss-of-function of Lssaw1 downregulates the expression of many adaxial genes but upregulates abaxial genes. LsSAW1 binds to the promoter region of the adaxial gene ASYMMETRIC LEAVES 1 (LsAS1) to upregulate its expression. Overexpression of LsAS1 compromised the effects of Lssaw1 on heading. LsSAW1 also binds to the promoter region of the abaxial gene YABBY 1 (LsYAB1), but downregulates its expression. Overexpression of LsYAB1 led to bending leaves in LsSAW1 genotypes. LsSAW1 directly interacts with KNOTTED 1 (LsKN1), which is necessary for leafy heads in lettuce. RNA-seq data showed that LsSAW1 and LsKN1 exert antagonistic effects on the expression of thousands of genes. LsSAW1 compromises the ability of LsKN1 to repress LsAS1. Our results suggest that downregulation or loss-of-function of adaxial genes and upregulation of abaxial genes allow for the development of leafy heads.

Up-regulation of LsKN1 promotes cytokinin and suppresses gibberellin biosynthesis to generate wavy leaves in lettuce.

Lettuce (Lactuca sativa) is one of the most popular vegetables worldwide, and diverse leaf shapes, including wavy leaves, are important commercial traits. In this study, we examined the genetics of wavy leaves using an F2 segregating population, and identified a major QTL controlling wavy leaves. The candidate region contained LsKN1, which has previously been shown to be indispensable for leafy heads in lettuce. Complementation tests and knockout experiments verified the function of LsKN1 in producing wavy leaves. The LsKN1∇ allele, which has the insertion of a transposon and has previously been shown to control leafy heads, promoted wavy leaves in our population. Transposition of the CACTA transposon from LsKN1 compromised its function for wavy leaves. High expression of LsKN1 up-regulated several key genes associated with cytokinin (CK) to increase the content in the leaves, whereas it down-regulated the expression of genes in the gibberellin (GA) biosynthesis pathway to decrease the content. Application of CK to leaves enhanced the wavy phenotype, while application of GA dramatically flattened the leaves. We conclude that the changes in CK and GA contents that result from high expression of LsKN1 switch determinate cells to indeterminate, and consequently leads to the development of wavy leaves.

LsNRL4 enhances photosynthesis and decreases leaf angles in lettuce.

Lettuce (Lactuca sativa) is one of the most important vegetables worldwide and an ideal plant for producing protein drugs. Both well-functioning chloroplasts that perform robust photosynthesis and small leaf angles that enable dense planting are essential for high yields. In this study, we used an F2 population derived from a cross between a lettuce cultivar with pale-green leaves and large leaf angles to a cultivar with dark-green leaves and small leaf angles to clone LsNRL4, which encodes an NPH3/RPT2-Like (NRL) protein. Unlike other NRL proteins in lettuce, the LsNRL4 lacks the BTB domain. Knockout mutants engineered using CRISPR/Cas9 and transgenic lines overexpressing LsNRL4 verified that LsNRL4 contributes to chloroplast development, photosynthesis and leaf angle. The LsNRL4 gene was not present in the parent with pale-green leaves and enlarged leaf angles. Loss of LsNRL4 results in the enlargement of chloroplasts, decreases in the amount of cellular space allocated to chloroplasts and defects in secondary cell wall biosynthesis in lamina joints. Overexpressing LsNRL4 significantly improved photosynthesis and decreased leaf angles. Indeed, the plant architecture of the overexpressing lines is ideal for dense planting. In summary, we identified a novel NRL gene that enhances photosynthesis and influences plant architecture. Our study provides new approaches for the breeding of lettuce that can be grown in dense planting in the open field or in modern plant factories. LsNRL4 homologues may also be used in other crops to increase photosynthesis and improve plant architecture.

The upregulated LsKN1 gene transforms pinnately to palmately lobed leaves through auxin, gibberellin, and leaf dorsiventrality pathways in lettuce.

Leaf shape represents a vital agronomic trait for leafy vegetables such as lettuce. Some lettuce cultivars produce lobed leaves, varying from pinnately to palmately lobed, but the genetic mechanisms remain unclear. In this study, we cloned one major quantitative trait locus (QTL) controlling palmately lobed leaves. The candidate gene, LsKN1, encodes a homeobox transcription factor, and has been shown previously to be critical for the development of leafy heads in lettuce. The LsKN1 allele that is upregulated by the insertion of a transposon promotes the development of palmately lobed leaves. We demonstrated that LsKN1 upregulated LsCUC2 and LsCUC3 through different mechanisms, and their upregulation was critical for the development of palmately lobed leaves. LsKN1 binds the promoter of LsPID to promote auxin biosynthesis, which positively contributes to the development of palmately lobed leaves. In contrast, LsKN1 suppresses GA biosynthesis to promote palmately lobed leaves. LsKN1 also binds to the promoter of LsAS1, a dorsiventrality gene, to downregulate its expression. Overexpression of the LsAS1 gene compromised the effects of the LsKN1 gene changing palmately to pinnately lobed leaves. Our study illustrated that the upregulated LsKN1 gene led to palmately lobed leaves in lettuce by integrating several downstream pathways, including auxin, gibberellin, and leaf dorsiventrality pathways.

Alternative splicing triggered by the insertion of a CACTA transposon attenuates LsGLK and leads to the development of pale-green leaves in lettuce.

Lettuce (Lactuca sativa) is one of the most important vegetable crops in the world. As a leafy vegetable, the polymorphism of lettuce leaves from dark to pale green is an important trait. However, the genetic and molecular mechanisms underlying such variations remain poorly understood. In this study, one major locus controlling the polymorphism of dark- and pale-green leaves in lettuce was identified using genome-wide association studies (GWAS). This locus was then fine mapped to an interval of 5375 bp on chromosome 4 using a segregating population containing 2480 progeny. Only one gene, homologous to the GLK genes in Arabidopsis and other plants, is present in the candidate region. A complementation test confirmed that the candidate gene, LsGLK, contributes to the variation of dark- and pale-green leaves. Sequence analysis showed that a CACTA transposon of 7434 bp was inserted 10 bp downstream of the stop codon of LsGLK, followed by a duplication of a 1826-bp fragment covering exons 3-6 of the LsGLK gene. The transposon insertion did not change the expression level of the LsGLK gene. However, because of alternative splicing, only 6% of the transcripts produced from the transposon insertion were wild-type transcripts, which led to the production of pale-green leaves. An evolutionary analysis revealed that the insertion of the CACTA transposon occurred in cultivated lettuce and might have been selected in particular cultivars to satisfy the diverse demands of consumers. In this study, we demonstrated that a transposon insertion near a gene may affect its splicing and consequently generate phenotypic variations.

Non-vernalization requirement in Chinese kale caused by loss of BoFLC and low expressions of its paralogs.

KEY MESSAGE: We identified the loss of BoFLC gene as the cause of non-vernalization requirement in B. oleracea. Our developed codominant marker of BoFLC gene can be used for breeding program of B. oleracea crops. Many species of the Brassicaceae family, including some Brassica crops, require vernalization to avoid pre-winter flowering. Vernalization is an unfavorable trait for Chinese kale (Brassica oleracea var. chinensis Lei), a stem vegetable, and therefore it has been lost during its domestication/breeding process. To reveal the genetics of vernalization variation, we constructed an F2 population through crossing a Chinese kale (a non-vernalization crop) with a kale (a vernalization crop). Using bulked segregant analysis (BSA) and RNA-seq, we identified one major quantitative trait locus (QTL) controlling vernalization and fine-mapped it to a region spanning 80 kb. Synteny analysis and PCR-based sequencing results revealed that compared to that of the kale parent, the candidate region of the Chinese kale parent lost a 9,325-bp fragment containing FLC homolog (BoFLC). In addition to the BoFLC gene, there are four other FLC homologs in the genome of B. oleracea, including Bo3g005470, Bo3g024250, Bo9g173370, and Bo9g173400. The qPCR analysis showed that the BoFLC had the highest expression among the five members of the FLC family. Considering the low expression levels of the four paralogs of BoFLC, we speculate that its paralogs cannot compensate the function of the lost BoFLC, therefore the presence/absence (PA) polymorphism of BoFLC determines the vernalization variation. Based on the PA polymorphism of BoFLC, we designed a codominant marker for the vernalization trait, which can be used for breeding programs of B. oleracea crops.

Two TOBAMOVIRUS MULTIPLICATION 2A homologs in tobacco control asymptomatic response to tobacco mosaic virus.

The most common response of a host to pathogens is arguably the asymptomatic response. However, the genetic and molecular mechanisms responsible for asymptomatic responses to pathogens are poorly understood. Here we report on the genetic cloning of two genes controlling the asymptomatic response to tobacco mosaic virus (TMV) in cultivated tobacco (Nicotiana tabacum). These two genes are homologous to tobamovirus multiplication 2A (TOM2A) from Arabidopsis, which was shown to be critical for the accumulation of TMV. Expression analysis indicates that the TOM2A genes might play fundamental roles in plant development or in responses to stresses. Consistent with this hypothesis, a null allele of the TOM2A ortholog in tomato (Solanum lycopersicum) led to the development of bent branches and a high tolerance to both TMV and tomato mosaic virus (ToMV). However, the TOM2A ortholog in Nicotiana glauca did not account for the asymptomatic response to TMV in N. glauca. We showed that TOM2A family is plant-specific and originated from Chlorophyte, and the biological functions of TOM2A orthologs to promote TMV accumulation are highly conserved in the plant kingdom-in both TMV host and nonhost species. In addition, we showed that the interaction between tobacco TOM1 and TOM2A orthologs in plant species is conserved, suggesting a conserved nature of TOM1-TOM2A module in promoting TMV multiplication in plants. The tradeoff between host development, the resistance of hosts to pathogens, and their influence on gene evolution are discussed. Our results shed light on mechanisms that contribute to asymptomatic responses to viruses in plants and provide approaches for developing TMV/ToMV-resistant crops.

Transcriptome regulation of carotenoids in five flesh-colored watermelons (Citrullus lanatus).

BACKGROUND: Fruit flesh color in watermelon (Citrullus lanatus) is a great index for evaluating the appearance quality and a key contributor influencing consumers' preferences. But the molecular mechanism of this intricate trait remains largely unknown. Here, the carotenoids and transcriptome dynamics during the fruit development of cultivated watermelon with five different flesh colors were analyzed. RESULTS: A total of 13 carotenoids and 16,781 differentially expressed genes (DEGs), including 1295 transcription factors (TFs), were detected in five watermelon genotypes during the fruit development. The comprehensive accumulation patterns of carotenoids were closely related to flesh color. A number of potential structural genes and transcription factors were found to be associated with the carotenoid biosynthesis pathway using comparative transcriptome analysis. The differentially expressed genes were divided into six subclusters and distributed in different GO terms and metabolic pathways. Furthermore, we performed weighted gene co-expression network analysis and predicted the hub genes in six main modules determining carotenoid contents. Cla018406 (a chaperone protein dnaJ-like protein) may be a candidate gene for ß-carotene accumulation and highly expressed in orange flesh-colored fruit. Cla007686 (a zinc finger CCCH domain-containing protein) was highly expressed in the red flesh-colored watermelon, maybe a key regulator of lycopene accumulation. Cla003760 (membrane protein) and Cla021635 (photosystem I reaction center subunit II) were predicted to be the hub genes and may play an essential role in yellow flesh formation. CONCLUSIONS: The composition and contents of carotenoids in five watermelon genotypes vary greatly. A series of candidate genes were revealed through combined analysis of metabolites and transcriptome. These results provide an important data resource for dissecting candidate genes and molecular basis governing flesh color formation in watermelon fruit.

Whole-genome resequencing of 445 Lactuca accessions reveals the domestication history of cultivated lettuce.

Lettuce (Lactuca sativa) is an important vegetable crop worldwide. Cultivated lettuce is believed to be domesticated from L. serriola; however, its origins and domestication history remain to be elucidated. Here, we sequenced a total of 445 Lactuca accessions, including major lettuce crop types and wild relative species, and generated a comprehensive map of lettuce genome variations. In-depth analyses of population structure and demography revealed that lettuce was first domesticated near the Caucasus, which was marked by loss of seed shattering. We also identified the genetic architecture of other domestication traits and wild introgressions in major resistance clusters in the lettuce genome. This study provides valuable genomic resources for crop breeding and sheds light on the domestication history of cultivated lettuce.

Comparative Metabolomic Profiling of Citrullus spp. Fruits Provides Evidence for Metabolomic Divergence during Domestication.

Watermelon (Citrullus lanatus) is one of the most nutritional fruits that is widely distributed in the whole world. The nutritional compositions are mainly influenced by the genotype and environment. However, the metabolomics of different domestication status and different flesh colors watermelon types is not fully understood. In this study, we reported an extensive assessment of metabolomic divergence in the fruit flesh among Citrullus sp. and within Citrullus sp. We demonstrate that metabolic profiling was significantly different between the wild and cultivated watermelons, the apigenin 6-C-glucoside, luteolin 6-C-glucoside, chrysoeriol C-hexoside, naringenin C-hexoside, C-pentosyl-chrysoeriol O-hexoside, and sucrose are the main divergent metabolites. Correlation analysis results revealed that flavonoids were present in one tight metabolite cluster. The main divergent metabolites in different flesh-colored cultivated watermelon fruits are p-coumaric acid, 2,3-dihydroflavone, catechin, N-(3-indolylacetyl)-l-alanine, 3,4-dihydroxycinnamic acid, and pelargonidin o-hexoside. A total of 431 differentially accumulated metabolites were identified from pairwise comparative analyses. C. lanatus edible-seed watermelon (cultivars) and C. mucosospermus (wild) have similar fruit metabolic profiles and phenotypic traits, indicating that edible-seed watermelon may be a relative of wild species and a relatively primitive differentiation type of cultivated watermelon. Our data provide extensive knowledge for metabolomics-based watermelon improvement of Citrullus fruits meet their enhanced nutritive properties or upgraded germplasm utility values.

Genome sequence and evolution of Betula platyphylla.

Betula L. (birch) is a pioneer hardwood tree species with ecological, economic, and evolutionary importance in the Northern Hemisphere. We sequenced the Betula platyphylla genome and assembled the sequences into 14 chromosomes. The Betula genome lacks evidence of recent whole-genome duplication and has the same paleoploidy level as Vitis vinifera and Prunus mume. Phylogenetic analysis of lignin pathway genes coupled with tissue-specific expression patterns provided clues for understanding the formation of higher ratios of syringyl to guaiacyl lignin observed in Betula species. Our transcriptome analysis of leaf tissues under a time-series cold stress experiment revealed the presence of the MEKK1-MKK2-MPK4 cascade and six additional mitogen-activated protein kinases that can be linked to a gene regulatory network involving many transcription factors and cold tolerance genes. Our genomic and transcriptome analyses provide insight into the structures, features, and evolution of the B. platyphylla genome. The chromosome-level genome and gene resources of B. platyphylla obtained in this study will facilitate the identification of important and essential genes governing important traits of trees and genetic improvement of B. platyphylla.

Upregulation of a KN1 homolog by transposon insertion promotes leafy head development in lettuce.

Leafy head is a unique type of plant architecture found in some vegetable crops, with leaves bending inward to form a compact head. The genetic and molecular mechanisms underlying leafy head in vegetables remain poorly understood. We genetically fine-mapped and cloned a major quantitative trait locus controlling heading in lettuce. The candidate gene (LsKN1) is a homolog of knotted 1 (KN1) from Zea mays Complementation and CRISPR/Cas9 knockout experiments confirmed the role of LsKN1 in heading. In heading lettuce, there is a CACTA-like transposon inserted into the first exon of LsKN1 (LsKN1▽). The transposon sequences act as a promoter rather than an enhancer and drive high expression of LsKN1▽. The enhanced expression of LsKN1▽ is necessary but not sufficient for heading in lettuce. Data from ChIP-sequencing, electrophoretic mobility shift assays, and dual luciferase assays indicate that the LsKN1▽ protein binds the promoter of LsAS1 and down-regulates its expression to alter leaf dorsoventrality. This study provides insight into plant leaf development and will be useful for studies on heading in other vegetable crops.

Dissection of the domestication-shaped genetic architecture of lettuce primary metabolism.

Lettuce (Lactuca sativa L.) is an important vegetable crop species worldwide. The primary metabolism of this species is essential for its growth, development and reproduction as well as providing a considerable direct source of energy and nutrition for humans. Here, through investigating 77 primary metabolites in 189 accessions including all major horticultural types and wild lettuce L. serriola we showed that the metabolites in L. serriola were different from those in cultivated lettuce. The findings were consistent with the demographic model of lettuce and supported a single domestication event for this species. Selection signals among these metabolic traits were detected. Specifically, galactinol, malate, quinate and threonate were significantly affected by the domestication process and cultivar differentiation of lettuce. Galactinol and raffinose might have been selected during stem lettuce cultivation as an adaption to the local environments in China. Furthermore, we identified 154 loci significantly associated with the level of 51 primary metabolites. Three genes (LG8749721, LG8763094 and LG5482522) responsible for the levels of galactinol, raffinose, quinate and chlorogenic acid were further dissected, which may have been the target of domestication and/or affected by local adaptation. Additionally, our findings strongly suggest that human selection resulted in reduced quinate and chlorogenic acid levels in cultivated lettuce. Our study thus provides beneficial genetic resources for lettuce quality improvement and sheds light on the domestication and evolution of this important leafy green.